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Abstract 

We notice that looking for v e — * v T at the same time as v e — * Vn oscillations 
could significantly help to reduce the errors in the leptonic CP-violating phase 
5 measurement. We show how the v e — ► ("golden") and v e — > v T ("silver") 
transitions observed at an OPERA-like 2 Kton lead-emulsion detector at L = 
732 Km, in combination with the v e — * v„ transitions observed at a 40 Kton 
magnetized iron detector with a baseline of L = 3000 Km, strongly reduce the so- 
called (#i3, 6) ambiguity. We also show how a moderate increase in the OPERA- 
like detector mass (4 Kton instead of 2 Kton) completely eliminates the clone 
regions even for small values of #13. 
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1 Introduction 



and solar 0-[14] neutrino data are strongly supporting 
and can be easily accommodated in a 



The present atmospheric |I[ 
the hypothesis of neutrino oscillations |15]-(T| 
three family mixing scenario. 

Let the Pontecorvo-Maki-Nakagawa-Sakata Upmns matrix be the leptonic analogue 
of the hadronic Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix in its most con- 
ventional parametrization |19|: 
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with the short-form notation s 



sin Oij , Cij 



cos 9ij. Oscillation experiments are 
sensitive to the two neutrino mass differences AmJ 2 , Am| 3 and to the four parameters 
in the mixing matrix: three angles and the Dirac CP-violating phase, 5. 

In particular, data on atmospheric neutrinos are interpreted as oscillations of muon 
neutrinos into neutrinos that are not v e 's, with a mass gap that we denote by Am^. 
The corresponding mixing angle is close to maximal, sin 2 2^23 > 0.8, and |Am| 3 | 



is m 



the range 1.8 to 4 x 10" d eV 2 |20 . 



The recent SNO results for solar neutrinos ||T2|-||14[ favour the LMA-MSW |2T 



solution of the solar neutrino deficit with v e oscillations into active (z^, v. 



neutrino 

states. The corresponding squared mass difference, that in this parametrization should 

is ~ 10~ 5 -10 -4 eV 2 . Comprehensive analyses of the solar 



be identified with Am 2 2 , 



neutrino data, however, do not exclude the LOW-MSW solution 



with Am 
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10 7 eV 2 . In both cases, the corresponding mixing angle (#12) is large (albeit not 
maximal) . 

Finally, the LSND data [26, 27] would indicate a — > v e oscillation with a third, 
very distinct, neutrino mass difference: Am^ SJVO ~ 0.3 — 6 eV 2 . The LSND evidence 
in favour of neutrino oscillation has not been confirmed by other experiments so far 
[ 28] ; the MiniBooNE experiment |H| will be able to do it in the near future [|30|. In the 
absence of an independent confirmation of the LSND evidence, we restrict ourselves 
to the three neutrino mixing scenario (the impact of a Neutrino Factory in the case of 
four neutrino mixing has been discussed in full detail in PT]-P3"| ). 

These oscillation signals will be confirmed in ongoing and planned atmospheric and 
solar neutrino experiments, as well as in long baseline ones, with the latter being free of 
model-dependent estimations of neutrino fluxes. There is a strong case for going further 
in the fundamental quest of the neutrino masses and mixing angles, as a necessary step 
to unravel the fundamental new scale (s) behind neutrino oscillations. In particular, 
it is possible that in ten years from now no information whatsoever will be at hand 
regarding the #13 angle (the key between the atmospheric and solar neutrino realms, for 
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which the present bound is sin 2 (26 , i 3 ) < 1 x 10 1 , [p^| ) and the leptonic CP violating 
phase 5. 

An experimental set-up with the ambitious goal of precision measurement of the 
whole three-neutrino mixing parameter space is under study. This experimental pro- 
gramme consists of the development of a "Neutrino Factory" (high-energy muons de- 
caying in the straight section of a storage ring and producing a very pure and intense 
neutrino beam, 



35, 36 



and of suitably optimized detectors located far away from the 
neutrino source. The effort to prepare such very long baseline neutrino experiments 
will require a time period covering this and the beginning of the following decade. It is 
therefore of interest to look for the optimal conceivable factory-detectors combination. 
One of its main goals would be the discovery of leptonic CP violation and, possibly, its 
study |3-|0 



Previous analyses [41 



}43f on the foreseeable outcome of experiments 
at a Neutrino Factory have shown that the determination of the two still unknown 
parameters in the three-neutrino mixing matrix, #13 and 5, will be possible (if the 
LMA-MSW solution of the solar neutrino deficit is confirmed). The most sensitive 
method to study these topics is to measure the transition probabilities involving v e 



and v ei in particular v e {y^) 



WW 



This is what is called the "golden measurement 



at the neutrino factory" . 
intense vJvp) beams. 



Such a facility is indeed unique in providing high energy and 

the transitions of inter- 



Since these beams contain no v^xy^, 
est can be measured by searching for "wrong-sign" muons: negative (positive) muons 
appearing in a massive detector with good muon charge identification capabilities |i0|. 

An incredible amount of work has been devoted to this topic in the last few years: 
we refer the interested reader to [44|-[54[ and to the refs. therein for an overview of 
the status-of-the-art in all its different aspects; we address to f^-pl and refs. therein 



for a comparison of the physics reach of a conventional (super)beam and of a Neutrino 
Factory; eventually, we point out that in |61[] the idea of a v e beam originating from 
/?-decay (the so-called "/3-beam"), was advanced: it appears that the physics reach of 
such a beam is complementary to that of a conventional superbeam 

it has been noticed that the probability P v 



In 53 



?i 3 , S, E v ) for neutrinos at 
a fixed energy (and for a given baseline) computed for a given theoretical input pair 
(6*i3, 5) defines a continuous equiprobability curve in the (#13, S) plane. Therefore, for 
a fixed energy, a continuum of solutions reproduce the input probability. A second 
equiprobability curve is defined in this plane by the probability for antineutrinos at 
the same energy and with the same input parameters, -Pp a ^p /3 (^i3, S, E v ). The two 
equiprobability curves have, quite generally, two intersection points: the first of them 
at (6\ 3 , 5), the second at a different point (#13, S). It is the intersection of equiprobability 
curves from the neutrinos and antineutrinos that resolves the continuum degeneracy 
of solutions in the (#13, 5) plane, restricting the allowed values for #13 and 5 to the two 
regions around (6*i 3 , 5) and (6*13,5). This second intersection, however, introduces an 
ambiguity in the measurement of the physical values of #13 and 6. Different proposals 
have been suggested to solve this ambiguity: in |63| the ambiguity is solved by fitting 
at two different baselines at the same time; another possibility is an increase in the 



3 



energy resolution of the detector, j|6], ^1|] ; see also |64j] . 



New degeneracies have later been noticed [pq] , resulting from our ignorance of the 
sign of the Am 2 atm squared mass difference (by the time the Neutrino Factory will be 
operational) and from the approximate [#23,v"/2 — #23] symmetry for the atmospheric 
angle. 

In the first part of the paper we describe how the (6*13, 8) ambiguity arises in v e — > 
oscillation due to the equiprobability curves in the (#13, 5) plane at fixed neutrino en- 
ergy. We then extend our analysis showing how the same phenomenon can be observed 
in a real experiment: equal- number-of-events (ENE) curves for any given neutrino en- 
ergy bin appears and their intersections in the (#13, 5) plane explain how and where do 
"clone" regions arise. Our analysis is afterwards compared with the results of simula- 
tions with a realistic magnetized iron detector (studied in [EH]). 

We then propose to reduce the continuum degeneracy and to resolve the last am- 
biguity using two baselines with detectors of different design. We notice that muons 
proceeding from r decay when r's are produced via a v e — > v T transition show a differ- 
ent (#13,5) correlation from those coming from v e — > z/ M (first considered in [flU]). By 
using a near lead-emulsion detector, capable of the r-decay vertex recognition, we can 
therefore use the complementarity of the information from v e — > v T and from v e — > 
to solve the (#13,5) ambiguity. We find that the combination of the near emulsion 
detector and of a massive magnetized iron detector at L = 3000 Km could indeed help 
to achieve a good resolution in the (#13, 5) plane. 

In this paper we restrict ourselves to the (#13, £) ambiguity, by fixing 823 = 45° 
and by choosing a given sign for Am 2 atm (in the hypothesis that more information on 
the three neutrino spectrum will be available by the time the Neutrino Factory will 
be operational). However, we found that the core of our results do not depend on 
the sign of Am^ m : indeed, we have been carrying on simulations with the opposite 
sign, with similar results (i.e., we still observe how "clone" regions disappear due to 
the two-detector types combination; notice, however, that the location of the "clone" 
regions and all the details of the simulation do depend on the sign of Am^ m ). We 
remind that, if # 13 is not extremely small, the combined measurement of v e —>■ and 
v e — > v T transitions could significantly help in solving the [#23, tt/2 — #23] ambiguity. 

In Section || we present our analysis of the v e — > equiprobability curves in 
the (#13,5) plane; in Section |3| we introduce the corresponding equal- number-of-events 
(ENE) curves; in Section |4] we present a similar analysis for the v e — > v T oscillation 
probability and study the impact of "silver" wrong-sign muon events; in Section |5| we 
show our results for the combination of a near (L = 732 Km) OPERA-like detector 
and of a L = 3000 Km magnetized iron detector; in Section |6] we eventually draw our 
conclusions. In Appendix [A] a perturbative expansion in A8 of the formulae of Sect. ^ 
is presented; in Appendix [B] we report some useful formulae for r CC-interaction and 
decay. 
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V 



equiprobability curves in the (6^13, <5) plane 



We consider the v e — > transition, first. This channel has been shown to be the 
optimal one to measure simultaneously #13 and S at the Neutrino Factory in the context 
of three-family mixing, through the appearance of "wrong-sign" muons in the detector, 
| 4~0fl . It therefore deserves the nickname of "golden channel" . 

Following eq. (1) of p3] we get for the transition probability at second order in 



perturbation theory in 6^3, A Q /A atm , A /v4 and A L (see also |67|-||69||), 

A a t m L 



^J(#i3, 6) = X ± sin 2 (2# 13 ) + Y ± cos(# 13 ) sin(2# 13 ) cos \±5 
where ± refers to neutrinos and antineutrinos, respectively, and 

sin 2 (^ 3 )(^) 2 sin 2 (^) , 
sin(2fl 12 ) sin(2£ 23 ) (^) (%^) sin ) sin (^f) 



+ Z. 





x± 
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Y ± 




Z 


with Z = Z + = 1 


Bzp = \A^f A atm \ 


nally, A atm 


= Am 



(2) 



cos 2 (^ 23 )sin 2 (2^ 12 ) 



sin 



2 (f) 



Z_. In these formulae, A = ^j2Gpn e (expressed in eV 2 /GeV) and 
(with =p referring to neutrinos and antineutrinos, respectively). Fi- 
J2E U and A = Am%/2E U . 

The parameters #i 3 and 8 are the physical parameters that must be reconstructed by 
fitting the experimental data with the theoretical formula for oscillations in matter. In 
what follows, the other parameters have been considered as fixed quantities, supposed 
to be known by the time when the Neutrino Factory will be operational. In particular, 

" 4 eV 2 [p^] -|^5|, corresponding 



in the solar sector we fixed #i 2 = 33° and Am = 1.0 x 10" 
to the LMA region of the solar neutrino problem (accordingly to the recent SNO results 
f0| 0); in the atmospheric sector, 6 23 = 45° and Am 2 tm = 2.9 x 10~ 3 eV 2 J20| (notice 
that for 823 = 45° the [6 l 23 ,7r/2 — 6> 23 ] ambiguity is absent). Finally, we considered 
a fixed value for the matter parameter, A = 1.1 x 10~ 4 eV 2 /GeV for L < 4000 km 
and A = 1.5 x 10~ 4 eV 2 /GeV for L > 4000 km, obtained by using the average matter 
density alongside the path for the chosen distance computed with the Preliminary 
Earth Model |70fl . For simplicity, we have not included errors on these parameters^. 

Eq. ([[]) leads to an equiprobability curve in the plane (#i 3 ,5) for neutrinos and 
antineutrinos of a given energy: 



(3) 



4 It has been shown in |33| that the inclusion of the foreseeable uncertainties on these parameters 
does not modify the results on the #13 and 6 measurements in a relevant manner. 
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We can solve eq. (|3|) for 5: 



1 2 J' y ± cos(0 13 )sin(20 13 ) ' [) 



It is useful to introduce the following functions: 



f(a a \ _ sin 2 (20i 3 ) -sin 2 (26>i 3 ) 
JI^ISj^IS; — cos6» 13 sin(26» 13 ) 



t ^1 3 A 3 ) 



cos 0i3 sin(20i3) 
cos 6»i3 sin(26»i 3 ) ' 



(5) 



with the obvious limit f(9 13 ,9 r3 ) = and g(9i 3 , 9 13 ) = 1. 
Eq. (|J) can then be written as: 

cos (±5 - ^f^j = R ± f(9 13 , 9 13 ) + cos (±5 - ^f^j g(9 13 , 9 13 ) . (6) 

Eq. (^) is particularly illuminating: it describes a family of two branches curves in 
the plane (9\ 3 , 5) for the neutrinos and a second family of two branches curves for 
the antineutrinos. The dependence on the neutrino energy resides in A atrn and in the 
ratio R± = X±/Y±, whereas the dependence on the angle is factorized in the two 
^-dependent functions / and g. 

It is helpful to introduce the parameter A#: 

#13 = #13 + A® ? 



constrained by the bound 



cos ( ±5 



AatmL 



< 1 



(7) 



The allowed region for A9 depend on the input parameters (9i 3 ,5), on the neutrino 
energy and on the baseline. In Fig. [I], we compute the allowed values of A9 as a 
function of 9 13 (9 13 G [0°, 13°]) at three different baselines, L = 732,3000 and 7332 
Km, for 6 = and E v = 38 GeV, by numerically solving eq. (0). 

It may be noticed that, for (almost) every value of #13 in the considered range, two 
different regions of allowed values for A# exist. The first region corresponds to A9 ~ 0, 
whereas the second corresponds to large negative values for A9. In this region, #13 is 
negative: this region is therefore unphysical when the sign of the mass differences and 
of the various angles are defined in an appropriate way |52 |. We concentrate hereafter 
on the tiny region around A9 ~ 0. 

In Fig. we present the equiprobability curves for P e + (the upper row) and P e ~ (the 
lower row) in the (A#, 5) plane, at L = 732, 3000 and 7332 Km for different values of 
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732 Km 



3000 Km 



7332 Km 





2 4 6 



Figure 1: Allowed region in A8 as a function of 9^ for E u = 38 GeV at three different 
distances: L = 732, 3000 and 7332 Km. 



the neutrino energy in the range E v G [5, 50] GeV. The input values are #13 = 5° and 
5 = 60°. In the upper row (neutrinos), it can be seen that all the equiprobability curves 
intersect in A8 = 0°,5 = 60° (namely, #13 = #13 and 5 — 5). However, notice that the 
equiprobability curve for a given neutrino energy intersects the curves corresponding 
to a different neutrino energy in a second point, at positive A9 and negative 5. This 
second intersection depends on the energies of the two curves. In the upper branch 
of the neutrino equiprobability curves, no second intersection is observed, for these 
particular values of the input parameters. The results of Fig. || may be understood 
with the help of a perturbative expansion of eqs. (|5]) and (|6|) in terms of powers of Ad 
(always possible in the allowed region, for #13 large enough). Details on this expansion 
can be found in App. [A[ 

We can draw some conclusion from what observed in Fig. |2| and from the previous 
considerations on the energy dependence of the equiprobability curves. In particular, 
it is to be expected that by fitting experimental data for neutrinos only it should be 
quite difficult to determine the physical parameters (0i3,5) with good accuracy. We 
expect, instead, that the fitting procedure will identify a low x 2 region whenever the 
family of equiprobability curves are not well separated, within the experimental energy 
resolution. In particular, at short distance (L = 732 Km) it is to be expected a good 
determination of #13 (notice that A6 is generally less than 2°) and no determination 
whatsoever of the CP-violating phase 5. At the intermediate distance, L = 3000 Km, 
the equiprobability curves for neutrinos (for this particular set of input parameters) do 
not depend strongly on the energy in the upper branch, whereas a larger separation 
can be seen in the lower branch. Therefore, we expect a low x 2 region alongside 
the upper branch of the equiprobability curves spanning from around the single point 
corresponding to the physical parameters (at A9 = 0, 5 = 5) to the (diluted) region 
where the curves show the second intersection (that by periodicity in the 5 axis happens 
to be in the lower branch). 
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1.2 2.4 



Figure 2: Equiprobability curves in the (A8, 5) plane, for 6* 13 = 5°, 5 = 60°, E v G [5, 50] 
GeV and L = 732, 3000 and 7332 Km. The upper row represents equiprobability curves 
for neutrinos, the lower row for antineutrinos. The dashed line is E u = 5 GeV, the 
solid line is E v = 45 GeV; the dotted lines lie in between these two. 



Finally, at large distance (L = 7332 Km) we expect a low x 2 region in the region 
AO < 0.5° and 5 ~ 5 ± 100° (notice that the small spread in the variable 9 for this 
baseline is in agreement with Fig. p]). 

A great improvement in the reconstruction of the physical parameters from the ex- 
perimental data is achievable using at the same time neutrino and antineutrino data. 
This can be seen in Fig. [|, where the equiprobability curves for neutrinos and antineu- 
trinos (for the same input parameters as in Fig. |2|) have been superimposed. At short 
distance, the two family of equiprobability curves overlap for any value of 5, and no 
improvement is to be expected. However, at the intermediate distance the equiproba- 
bility curves for neutrinos and antineutrinos overlap only in the vicinity of the physical 
point A6 = 0, 5 = 5 and in the region of the second intersection, whereas in the in- 
termediate region they are quite well separated, both in the upper and lower branch. 
We expect, in this case, that the fitting procedure of the whole set of neutrino and 
antineutrino data will identify two separate low \ 2 region, around the physical point 
and around the region where all the curves show the second intersection. This second 
allowed region in the parameter space was first observed in |63|] and subsequently con- 
firmed in [[|6|, [5(| Finally, at large distance we expect no significant improvement 



with respect to the previous case. 

Notice that these considerations can be drawn by looking at the equiprobability 
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curves for neutrinos and antineutrinos, only. We will see in the following section how 
the theoretical expectation is indeed reproduced in the "experimental data" . 



3 Number of "wrong-sign" muons in the detector 

The experimental information is not the transition probability but the number of 
muons with charge opposite to that of the muons circulating in the storage ring, that 
in the following will be often called "golden" muons. The events are then grouped in 
bins of energy, with the size of the energy bin depending on the energy resolution AE 
of the considered detector. In general, 

E * +AE , s +/ - « d$ v {v )(E) 

ATP „ . . TP\ TD± TP O . _ X\ Ve(Ve)\ I o 



N^(d 13 ,5) = / dEa MV)i) (E)P±(E,e 13 ,S) 



dE 



is the number of wrong-sign muons in the i-th energy bin for the input pair (# 13 , 5); E 
is the neutrino (antineutrino) energy^. The charged current neutrino and antineutrino 
interaction rates can be computed using the approximate expressions for the neutrino- 
nucleon cross sections on an isoscalar target, 

a uN ~ 0.67 x 1(T 42 x J^- x m 2 ; a m ~ 0.34 x 10" 42 x x m 2 . (9) 

GeV GeV 

In the laboratory frame the neutrino fluxes, boosted along the muon momentum vector, 
are given by: 

fir - J^f^ 1 -^) {K -2^(1-^03,)] 

T^K" 2B> (1-/3 cos <p)]} . (10) 



5 The neutrino energy can be reconstructed if the considered detector has a hadronic calorimeter 
capable to measure the energy of the hadronic shower (Ehadr) in the vN CC interactions with good 
precision. 
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Here, (3 — J 1 — (m^/E^) 2 , is the parent muon energy, y = E v jE^ n M is the number 
of useful muons per year obtained from the storage ring and ip is the angle between the 
beam axis and the direction pointing towards the detector. In what follows, the fluxes 
have been integrated in the forward direction with an angular divergence (taken to be 
constant) 5(f ~ 0.1 mr. The effects of the beam divergence and the QED one-loop 
radiative corrections to the neutrino fluxes have been properly taken into account in 
7j| . The overall correction to the neutrino flux has been shown to be of 0(0.1%). 

In the same approximations as for eq. ([[]), we get for the number of events per bin: 



Nf. 



It sin 2 (20i3 ) + Iyc cos 5 + I Y s sin 6 cos 13 sin(20i 3 ) + P z 



F x _ sin 2 (20 13 ) + 



Iyc cos S — Ivs sin 5 



cos0 13 sin(20 13 ) + 4, 



where we introduced a short-form notation for the following integrals: 



p 


= f^UEa MD /^f/ E) X ± (E), 




p 

1 vc 


= J^ E dEa Mp /^f/ E) Y ± (E) cos 


( A atm L \ 
\ 2 ) 


p 

r ± 


= J^ E dEa v ^f/ E) Y ± (E)sm 


( AatmL \ 

\ 2 ; 


p 


= !l;^ E dEa^/^^Z(E). 
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For a fixed energy bin and fixed input parameters (0i 3 , 5), we can draw a continuous 
curve of equal number of events in the (A0, 5) plane, 



iv; ± (0 13 ,5) = iv; ± (0 13 ,5), 

as it was the case for the transition probability, eq. (|3|). 
We therefore get an implicit equation in 5, 

F(5) = G(0 13 ,0 13 ,5), 

where 



(13) 



F(S) 

£(013,013,5) 



cos 5 ± yt^- sin 5 



(14) 



(15) 



7^ )f(0i3,e 13 ) + F(S)g(6 13 ,e l3 ) 



and /(0i 3 ,0i 3 ) and g(0i 3 ,0i 3 ) are the 0-dependent functions introduced in eq. (^). 
Solving for 5, 

5 = F- l [G(9 13 ,9 13 ,5)} (16) 
we get equal-number-of-events curves (ENE) in the ( A0, 5) plane, see Fig. f|. 
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Figure 4: Equal-number- of-events (ENE) curves in the (AO, 5) plane for neutrinos, for 
13 = 5°,5 = 60°, E v G [5,50] GeV and L = 732, 3000 and 7332 Km. The dashed 
line represents the first energy bin, E v G [0, 10] GeV, the solid line the last energy bin, 
E v G [40,50] GeV; the dotted lines lie in between these two. 



The problem arises in the reconstruction of the physical parameters from a data 
set consisting of some given number of events per bin, for a given number of bins 
(depending on the specific detector energy resolution). As it was the case for the 
equiprobability curves in the previous section, all the ENE curves intersect in the 
physical point (A9 = 0,5 = 5) and any given couple of curves intersect in a second 
point in the same region as in Fig. |2|. As it can be seen in Fig. [|, the second intersection 
differs when considering different couples of curves, but lies always in a restricted area of 
the (A9, 5) plane, the specific location of this region depending on the input parameters 
(9i3,5), see eqs. ( p5|) and (p6|). The x 2 analysis of the data will therefore identify two 
allowed regions: the "physical" one (around the physical value, 9i3,5) and the "clone" 
solution, spanning all the area where a second intersection between any two ENE curves 
occurs. This is the source of the ambiguity pointed out in [B3T]. 



In the remaining of this section, we apply the analysis in energy bins of [fH], |63 
Let iVA, be the total number of wrong-sign muons detected when the factory is run in 
polarity p = grouped in energy bins specified by the index i, and three possible 

distances, A = 1, 2, 3 (corresponding to L = 732 Km, L = 3000 Km and L = 7332 Km, 
respectively). In order to simulate a typical experimental situation we generate a set 
of "data" n x p as follows: for a given value of the oscillation parameters, the expected 
number of events, N x p , is computed; taking into account backgrounds and detection 
efficiencies per bin, b x p and e x p , we then perform a gaussian (or poissonian for N x p < 10 
events) smearing to mimic the statistical uncertainty: 



A 



i:p 



Smear(A^ p + tf ) . (17) 



Finally, "data" are fitted to the theoretical expectation as a function of the neutrino 
parameters under study, using a gaussian x 2 minimization: 



x - K x 
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where 5nf p is the statistical error for nf p (errors on background and efficiencies are 
neglected) or a poissonian \ 2 minimization: 



-2EE 




(19) 



whenever events are Poisson-distributed around the theoretical values (see |60| and 
refs. therein). We verified that the fitting of theoretical numbers to the smeared 
("experimental") ones is able to reproduce the values of the input parameters (the 
best fit always lies within a restricted region around 9\3 , 8) . 

The following "reference set-up" has been considered: neutrino beams resulting 
from the decay of 2 x 10 20 /i +, s and /i~'s per year in a straight section of an = 50 
GeV muon accumulator. An experiment with a realistic 40 Kton detector of magnetized 
iron and five years of data taking for each polarity is envisaged. Detailed estimates 
of the corresponding expected backgrounds and efficiencies have been included in the 
analysis, following [J66| . Notice that this set-up is exactly the same of [fE], j33|j . 

In the first row of Figs. [||], we present the results of the fit to five bins of data 
for decaying muons of one single polarity, fi~. The energy resolution of the detector is 
AE U = 10 GeV. In all cases we observe the pattern depicted in the previous section, 
with a good determination of #13 and an extremely poor determination of 5. In the 
second row we fit to five bins of data for decaying muons of both polarities. The 
results follow again the theoretical analysis of this and of the previous section and 
are in perfect agreement with what presented in [63]. In particular, notice how at the 



intermediate distance it is possible now to reconstruct 5 with an error of tens of degrees 
in two separate regions of the parameter space. 
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Figure 5: #5.5, 50 and 55 % C.L. contours resulting from a \ 2 fit of 9\s and 5, for 
9 is = 1° and 5 = 0° (up) and 5 = 90° (down). For each considered input parameters 
couple, the two rows represent: from left to right, L = 732, 3000 and 7332 Km; the 
upper row is iV* + , i = 1, . . . , 5; the lower row is iV* + and iV* ,i — 1, ... ,5. 
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Figure 6: 68.5, 90 and 99 % C.L. contours resulting from a \ 2 fit of 9^ and 5, for 
# 13 = 5° and 5 = 0° (up) and 5 = 90° (down). For each considered input parameters 
couple, the two rows represent: from left to right, L = 732, 3000 and 7332 Km; the 



upper row is N l + } i 



5: the lower row is N* 



and N*,i 
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4 The v e — > v T equiprobability and ENE curves 

We present in this section the possibility to use a different channel, namely the v e — > v T 
oscillation probability, to improve the reconstruction of the physical parameters (#i 3 , 5) 
in combination with the results for the u e — > transition described in j|0], and in 
the previous section. 

The v e — > v T oscillation probability at second order in perturbation theory in # 13 , 
A Q /A atm , A Q /A and A L is: 



P±(0i3, 5) = XI sin^(2# 13 ) - Y± cos(# 13 ) sin(2fl 13 ) cos ±5 



A 



atm J 



where ± refers to neutrinos and antineutrinos, respectively, and 



(20) 



Y ± 



sin 



m ■ 

sin(20 12 ) sin(20 23 ) (^) sin (*f) sin (i 



sin 2 (9 23 )sin 2 (2# 12 ) 



sm 



2 / al\ 



V 2 y ' 



(21) 



ZI. Notice that X£ and Z r differs from the corresponding coefficients 



with Z T = Z\ 

for the v e — > transition for the cos 823 *-* sin 623 exchange, only. The Y± term is 
identical for the two channels, but it appears with an opposite sign. This sign difference 
in the F-term is crucial, as it determines a different shape in the (AO, 5) plane for the 
two sets of equiprobability curves. 
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Figure 7: Equiprobability curves in the (AO, 5) plane, for O13 = 5°, 5 = 60° , E v G [5, 50] 
GeV and L = 732 Km for the u e and v e — > v T oscillation (neutrinos on the left, 

antineutrinos on the right). 
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In Fig. [7], we superimposed the equiprobability curves for the v e — > v T and v e — > 
oscillations at a fixed distance, L = 732 Km, with input parameters 9^ = 5° and 
5 = 60°, for different values of the energy, E u e [5, 50] GeV. The effect of the different 
sign in front of the Y-term in eqs. (0) and (|20|) can be seen in the opposite shape in 
the (#i3, S) plane of the v e — > v T curves with respect to the v e — > ones. Notice that 
all the curves of both families met in the "physical" point, #13 = #13, 5 = 5, and that 
now three would-be "clone" regions (i.e., the spread regions where the intersections of 
any given couple of equiprobability curves lie) can be seen. 

As a final comment we signal that, if 6 , 13 is not extremely small (in such a way 
that the X± and X± terms dominate over the Y± terms in eqs. (|TJ) and fl2"0p), the 
combined measurement of v e —>■ and v e — > v T transitions could in principle solve the 
[6> 23 , 7r/2 - 23 ] ambiguity. 

To follow the line of reasoning adopted for the v e — > channel, we should now 
discuss how the v e — > v T channel can be used in a realistic experiment. A number of 
modifications with respect to the case of the "golden" channel should be taken into 
account. 

First, the approximate expressions for the neutrino- nucleon cross section on an 
isoscalar target, eq. @, are no longer appropriate in the case of a v T CC interaction 
inside the detector. In this case we used the reported values for the v T N cross-section 
|72fl that have been applied in the CHORUS and OPERA experiment to compute the 



expected number of CC r events. The considered cross-section includes r mass effects 
in the DIS region following |73[] (see App. || for details), as well as the elastic and 
quasi-elastic contributions to the cross section. 

Second, the r will decay in flight into a muon of the same charge and two neutrinos, 



with a branching ratio BR(t — > fi) ~ 0.17, |19]. This "silver" wrong-sign muon is the 
experimental signal we are looking for, to be identified and to be separated from the 
"golden" wrong-sign muons originated from CC interactions^. The first tool to 
distinguish the two sets of wrong-sign muons is their different energy distribution (see 
App. [B] for details on the differential decay rate). It has been shown in |>6| that in 
the magnetized iron detector considered in the previous section, muons from r decay 
cannot be distinguished from the main background represented by muons from charmed 
mesons decay by means of kinematical cuts. In order to take advantage of this channel, 
we should therefore use a different kind of detector: for this reason we concentrate in 
the remaining of the paper on a lead-emulsion detector, where the observation of the 
r decay vertex allows to distinguish "golden" and "silver" wrong-sign muons, and 
the latter from the charmed mesons decay background. We must mention that the 



v e — > v T oscillations were previously considered in for a liquid argon detector, 
using kinematical cuts to identify r's. It could be of interest to explore further the 
possibility of using "silver" muon events in such a detector to reduce or eliminate the 



6 We adopt the nick-name of "silver" muon events due to the lesser statistical significance with 
respect to "golden" ones. 
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(6*13,5) ambiguity. 

In what follows, we consider an OPERA-like detector with a mass of 2 Kton and 
spectrometers capable of muon charge identification (see the OPERA proposal for 
details, [|74[]) located at L = 732 Km down the neutrino source (obviously, in the back 



of our mind we are thinking of the CNGS set-up). The results for both the "golden" 
and the "silver" channel at the near emulsion detector will be combined with results 
for the "golden" channel obtained with the magnetized iron detector located at the 
optimal distance for the measurement of leptonic CP violation, L = 3000 Km. 

In this paper, we will first restrict ourselves to an ideal OPERA-like detector with 
perfect efficiency and no background. Afterwards, we take into account the realistic 
estimates of the energy-dependent reconstruction efficiency and of the most relevant 
backgroundsQ Eventually, we will consider how an increase in the detector mass or an 
improvement on the signal/noise ratio affects our results. 

Schematically, starting from a positive charged muon in the storage ring, "silver" 
muons are obtained by the following chain: 



v e -> v T —> t~ -> \L 



whereas "golden" muons come from: 




If we group the events in bins of the final muon energy E^, with the size of the energy 
bin depending on the energy resolution AE^ of the considered detector, the number of 
"golden" muons in the i-th energy bin for the input pair (#13, 5) and for a parent muon 
energy E^ is: 

KM 3 j) = { j^agggj » p±(km » diM 'f E E ; ,B, ' ) } E ^ AE " (22) 

(remember that P ± is the oscillation probability for neutrinos and antineutrinos, re- 
spectively, see Sect. 0), whereas the number of "silver" muons in the i-th energy bin 

7 A dedicated careful analysis of the "silver" muons reconstruction efficiency and of the background 
as a function of the neutrino energy for this specific detector is currently under progress, [j7{|. A key 
issue is the maximum affordable amount of charge discrimination, both for the emulsion and the 
magnetized iron detector. Estimates can be found in J6^ | for the magnetized iron detector and in 
Fig. 86 of [Q for the emulsion detector. 
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Figure 8: Comparison of the number of "golden" muons (lightest bars) and "silver" 
muons (darkest bars), for a parent muon energy E^ = 50 GeV and input parameters 
#13 = 5°, 5 = 60° and a near detector (L = 732 Km) with a mass of 2 Kton and perfect 
reconstruction efficiency for both channels. 



is: 



dN^(E^,E T ) da UT{PT) (E T , E v ) 



dE„ 



dE T 



p±(E u ,e 13 ,5) ® 



d$> Ve (p e ){E v , E^ 



dE tJ 



(23) 



In both equations, <8> stands for a convolution integral on the intermediate energy: for 
example, 



BR{t -> n) 



dN^(E^E T ) da Ur{VT) {E T) E u 



dE„ 



dE T 



gives the number of muons in the i-th bin in the final muon energy E^, for a given 
neutrino energy E v (see App. ^| for details). In Fig. |8| we present a direct comparison 
of "golden" and "silver" muons grouped in five energy bins with AE^ = 10 GeV, 
for a parent muon energy = 50 GeV with input parameters #13 = 5°, 5 = 60°. 
We consider here a near detector (L = 732 Km) with a mass of 2 Kton and perfect 
reconstruction efficiency for both channels. Notice that in Fig. |S| we have not included 
quasi-elastic and resonance contributions to the uN cross-section that could enhance 
r production at low neutrino energy. 

The total number of events for these parameters are iVf_ ~ 700 and N* ~ 30. The 
strong reduction in the number of "silver" muons with respect to the "golden" muons 
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Figure 9: Equal-number- of- event-curves in the (AO, 5) plane, for E v G [5,50] GeY and 
and L = 732 Km in the case of v e — > v T oscillation (neutrinos on the left, antineutrinos 
on the right), for # 13 = 5° and 5 = 60°. 



Following the same procedure used to get the "golden" muons ENE curves presented 
in Fig. |] we can compute ENE curves for "silver" muons. These curves are reported 
in Fig. ^ in the case of L = 732 Km. 



In Fig. [10] we superimpose "golden" and "silver" ENE curves for L = 732 Km and 
6*i3 = 5°, 5 = 60°. Notice how, as it was expected from the equiprobability curves 
analysis, the two sets of curves have opposite concavity in the ( AO, 5) plane. As for the 
equiprobability curves, all lines met in the "physical" point. Therefore, a combined x 2 
analysis of the two sets of data should present a well defined global minimum around 
the "physical" region, whereas the local minima situated in the three "clone" regions 
are considerably raised with respect to what presented in the previous section where 
only the "golden" muon signal was considered. 



5 Combination of "golden" and "silver" muon events 



We follow the analysis in energy bins outlined in the previous section and in p0| , |63| : 



we produce a theoretical data set (A^ ± , iV*±) for fixed input parameters #i 3 , S and then 
we smear it as in eq. ( |ITD to obtain an "experimental" data set (n 9 ±,n 8 ±). Finally, 
"experimental" data are fitted and 68.5, 90 and 99 % C.L. contours in the (#13, 5) plane 
are drawn. 

In Fig. |1J] we present the results of this analysis comparing two different possibilities: 
in the upper row we combine two realistic magnetized iron detectors at L — 732 and 
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L = 3000 Km; in the lower row we combine an ideal OPERA-like detector at L = 732 
Km and a realistic magnetized iron detector at L — 3000 Km. 

First, we present our results for the combination of the two iron detectors (where 
only "golden" muons, N^±, can be used). We draw in each figure the contours for 
different input parameters: three values for #13 = 1°,6° and 11° and three values for 
the phase 5 = —90°, 0° and 90°. In each figure, therefore, fits to nine input parameter 
pairs (^13, 5) are shown: this has to be compared with Figs. |5| and |6] where in each plot 
the results of a fit to one single input parameter pair was presented. 
On the left, only five years of data taking for /i + circulating in the storage ring are 
considered. Notice that for any given input pair # 13 is always reconstructed within 
a 2° error; on the contrary, roughly any value for the CP-violating phase is allowed. 
The situation is drastically improved on the right, where five years of data taking for 
each muon polarization are considered. In this case, the phase 5 is reconstructed with 
a precision of tens of degrees for all values of the input parameters. Notice, however, 
how some "clone" region is still present at 90 % C.L. (e.g., for #13 = 1°,5 = 90° the 
small region around #13 = 2°, 5 = 150°; for #13 = 1°,5 = 0° the small region around 
#13 = 2.5°, 5 = —150°; see also Fig. [5]). These results can be easily understood in terms 
of the theoretical analysis of the equiprobability and equal-number-of-events curves for 
the v e — > channel of Sects. ^] and |3|. 

We present now results for the combination of a near emulsion detector (with 
both "golden" and "silver" muons, N 9 ± and N^±) and a not-so-far iron detector (with 
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Figure 11: 68.5, 90 and 99 % C.L. contours resulting from a \ 2 fit o/ #13 and 5, for 
#i3 = 1°,6° and 11°, and 5 = — 90°, 0° and90°,for the combination of two iron detectors 
(upper row) or one iron detector and one emulsion detector (lower row). Upper row: 
iron detector at L = 732 and 3000 Km: a) N 9 + ; b) N 9 + and N 9 ; Lower row: iron 
detector at L = 3000 Km and emulsion detector at L = 732 Km: a) N 9 + + N^+; b) 
N 9 +NUandN 9 . 



"golden" muons, only). On the left, again only five years of data taking for /i + cir- 
culating in the storage ring are considered. Notice that a significant reduction in the 
reconstruction errors on the phase 5 is already achieved. On the right, we simply add to 
the first five years of data taking for the /i + polarity further five years for the opposite 
polarity in the iron detector, only. We have not included a further five year operational 
time for the emulsion detector to take into account the mass decrease due to the brick 
removal in the first five year periodP]. Notice, however, that a quite relevant improve- 
ment with respect to the one-polarity two-detector types case (lower left) is achieved 
in the 5 reconstruction error. More important, an improvement with respect to the 
two-polarities one-detector type (upper right) can also be observed. In particular, the 
"clone" regions have completely disappeared, due to the combination of "silver" and 
"golden" muons with a different (^-dependent oscillation probability. The effect of the 

8 This is a quite conservative assumption, being the expected number of bricks to be removed 
looking for "silver" muons smaller than in the case of — > v T oscillations considered in the OPERA 
proposal, Q. 
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inclusion of "silver" muons can be seen in Fig. O, where we present the two-polarities 



two-detector types combination (Fig. [11], lower right) compared with the same com- 
bination but with only a "golden" muon signal. In this figure we can clearly see that 
"clone" regions are still present and that the near emulsion detector is too small to 
compete with a 40 Kton iron detector located at the same distance down the neutrino 
source. 
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Figure 12: 68.5, 90 and 99 % C.L. contours resulting from a \ 2 fit o/ #13 and 5, for 
813 = 1°, 6° and 11° , and 5 = —90°, 0° and 90°, for the combination of an iron detector 
at L = 3000 Km and an emulsion detector at L = 732 Km: (left) both "silver" and 
"golden" muon events are taken into account (same figure as in Fig. [77], lower right); 
(right) only "golden" muon events are taken into account. Five years of data taking 
for both polarities in the distant detector and only five years in the fi + polarity in the 
near detector have been considered. 



We also studied the effect of a change in the longest baseline, moving the iron 
detector from L = 3000 Km to L — 2000 Km, with no significant improvement. This 
was to be expected, since the former distance, L = 3000 Km has been found to be the 
optimal one to measure the CP- violating phase for the considered set-up, [P0|-pD[. 



We remind here that in [^3| it was shown that the optimal combination of two 
magnetized iron detectors is achieved for L = 3000 Km and L = 7332 Km: in that 
case, no "clone" region is present (see Fig. 6 of |63[). We prefer however to compare 
our results for the OPERA-like near detector and magnetized iron not-so-far detector 
option with the two magnetized iron detectors combination at the same distances. 

Up to this moment, we considered an ideal lead-emulsion detector with identification 
efficiency equal to 1 and no backgrounds. We should consider instead a real detector 
with a refined analysis of efficiencies and background for the two "wrong-sign muons" 
signals, u e — > and v e — > v T . A detailed simulation of the expected performance at an 
OPERA-like detector is under way ||75|| . However, a preliminary study has been done 



following the outline of the OPERA proposal |74| and of the recent progress report 
We first consider "golden" muons. In this case, it seems reasonable to consider 
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a 90 % average efficiency in the neutrino energy range E v G [5, 50] GeV. The most 
relevant background to this channel are "right-sign" muons with a wrong charge as- 
signment: with no improvement with respect to the present OPERA proposal, a level 
of 10~ 3 of charge misidentification is achieved (a level of 10 -6 was envisaged for the 
magnetized iron detector, see 
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On the other side, the dominant background at 
the magnetized iron detector coming from charmed mesons decay with no "right-sign" 
muon identification is here under control by looking for the decay vertex. 

Regarding "silver" muons, we have an unavoidable important reduction in the aver- 
age efficiency, mainly due to the detector design: the r decay vertex cannot be observed 
if it occurs inside the lead plates. An average efficiency for "silver" muons of 25 % only 
is achieved in the neutrino energy range E u e [5,50] GeV, [74]. In the case of "sil- 




Figure 13: Charmed mesons decay background and "silver" muons signal (shaded) for 
6 is = 11° and 5 = 90° in the OPERA-like emulsion detector as a function of the 
reconstructed neutrino energy (for 1 GeV energy bins). 



ver" muon events, the dominant background is the standard charmed mesons decay 
with no "right-sign" muon identification. Sub-dominant backgrounds coming from tt 
and K mesons decay can be kept at a reasonable level by imposing kinematical cuts 
as described in ||74fl . As an illustrative example, in Fig. [13] we present the charmed 
mesons decay background and the "silver" muon signal for #13 = 11° and 5 = 90° as 
a function of the neutrino energy (for 1 GeV energy bins). Notice how the "silver" 
muons signal peaks at lower energy with respect to the background, thus allowing us 
to use the u e —>■ v T oscillation to improve the (#13, 5) reconstruction. 



Eventually, we present in Fig. [TJ] the results for the combination of an iron detector 
at L — 3000 Km and of an emulsion detector at L = 732 Km with reasonable estimates 
for the reconstruction efficiencies for "golden" and "silver" muons and for the corre- 
sponding backgrounds. We can see that the results obtained with the two-detector 
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types combination strongly resembles those obtained with the reference two iron de- 
tectors set-up (Fig. [11], upper right). In particular, we notice a general improvement 
at the 68 % CL. but we still observe "clone" regions for small values of #13. 
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Figure 14: 68.5, 90 and 99 % CL. contours resulting from a \ 2 fit 0/ #13 and 5, for 
#13 = 1°, 6° and 11° , and 5 = —90°, 0° and 90°, for the combination of an iron detector 
at L = 3000 Km (with "golden" muon events) and an emulsion detector at L = 732 
Km (with "golden" and "silver" muon events). Realistic estimates for efficiencies and 
backgrounds at both the iron and the emulsion detector have been taken into account. 
Five years of data taking for both polarities in the distant detector and only five years 
in the /i + polarity in the near detector have been considered. 



Although the results of Fig. [L4] do not indicate a strong improvement with respect 
to the reference two iron detectors set-up, two comments are in order: first, we have 
been considering up to now an OPERA-like detector with the characteristics of the 
detector described in \74\ which is currently under construction and should be opera- 
tional by year 2006. We therefore study the effects on the physics reach of foreseeable 
improvements in the following. Secondly, it seems to us that the combination of two 
different detector types with different characteristics and systematics could be of the 
utmost importance in the project of a facility for precision measurements of the leptonic 
mixing matrix such as the Neutrino Factory. 

In Tab. [I] we schematically present which are the main advantages of an increase 
in the emulsion detector mass, or of a reduction of the main backgrounds, considering 
the presence or absence of "clone" regions and the maximal error in 5 that can be 
achieved. The main background for "silver" muon events B s u ver (originating from 
charmed mesons decay) is normalized to the value that can be deduced on the basis of 
the OPERA proposal. In this case we do not consider that significant improvements 
can be achieved and we restrict ourselves to a 10 to 20 % reduction. We believe that, 
on the other hand, the main background for "golden" muon events B go i den (originating 
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Table 1: Effects of a mass increase or of a "silver" (B s u ver ) or "golden" (B go id en ) main 
background reduction for the emulsion detector, for 6>i 3 = 1° and 5 = 90°. Whenever a 
"clone" region is present, no overall error on 5 is given. 

from charge misidentification) could be significantly reduced: we therefore present 
results for a IO -3 (the predicted level in the OPERA proposal) and for a 10~ 4 charge 
misidentification. Finally, having in mind that in the last years the scanning power 
increased by about a factor 10 every two years, we consider as realistic the assumption 
that, by the time the Neutrino factory will become operational, we will be able to 
scan a larger number of events, say a factor two. Therefore, we studied the physical 
reach of the Neutrino factory-detectors set-up as a function of the lead-emulsion mass, 
considering a linear increase in the mass from 2 Kton to 4 Kton. 

The outcome of this analysis is the following: 

• "Clone" regions disappear increasing the emulsion mass. For a 4 Kton detector 
no ambiguity is present. An improvement in the charge identification (i.e., a 
reduction of "golden" muons main background) implies generally a decrease in 
the statistical significance of "clone" regions. 

• When no "clone" regions are observed, a reduction of "golden" muons main 
background induces a shift in the 99 % C.L. allowed region for 5 towards higher 
values. In the optimal case we observe a gaussian-distributed region with A5 = 
±20°. 
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A 10 to 20 % reduction of the "silver" muons main background (charmed mesons 
decay) does not seem to improve significantly the previous results. 



We remind that this table should be interpreted as an indication of which kind of 
improvement in the emulsion detector helps more to improve the physics reach of the 
envisaged two detector types combination. 



6 Conclusions 



It was previously shown EH] that the "golden" wrong-sign muon signal at the Neutrino 
Factory can be extremely useful to measure simultaneously two of the parameters of 
the PMNS leptonic mixing matrix, #13 and 5. In |)3| was first noticed that degenerate 
regions in the (#13,5) parameter space occur in many cases, severely reducing the 
Neutrino Factory sensitivity to the CP-violating phase S. 

In the first part of the paper we describe how the (#13, 6) ambiguity arises in v e — > 
oscillation due to the equiprobability curves in the (#13, 5) plane at fixed neutrino en- 
ergy. We then extend our analysis showing how the same phenomenon can be observed 
in a real experiment: equal- number-of-events (ENE) curves for any given neutrino en- 
ergy bin appears and their intersections in the (#13, 5) plane explain how and where do 
"clone" regions arise. Our analysis is subsequently compared with the results of simu- 
lations with a realistic magnetized iron detector, already studied in previous papers. 

In the second part of the paper, we propose the use of v e — > v T oscillation to 
solve the (9 13 , 6) ambiguity problem. We apply the arguments presented in the case of 
v e v p, and we show how v e — > v T actually give a complementary information that can 
be used to improve the measurement of the phase S. However, to take full advantage of 
this "silver" channel we should use a detector of a different kind to distinguish muons 
originated from r decay from the so-called "golden" muons. In the last part of the 
paper we therefore perform simulations using the combination of a realistic magnetized 
iron detector at L = 3000 Km and of an OPERA-like lead-emulsion detector located 
at L — 732 Km from the Neutrino Factory. The r decay vertex recognition in the 
emulsion detector allow us to separate "silver" from "golden" muon events and to 
strongly reduce the effect of the charmed mesons decay background. 

We first use an ideal lead-emulsion detector with the mass currently envisaged by the 
OPERA Collaboration (2 Kton), perfect reconstruction efficiency and no background, 
showing how the "clone" regions disappear for any value of #13 > 1°. Eventually, 
we consider a realistic estimate of the reconstruction efficiency and of the main back- 
grounds both for "golden" and "silver" muon events at the emulsion detector, following 
the outline of the OPERA proposal |[74|| . In this case, "clone" regions for #13 ~ 1° do 
appear and our results are similar to those that we obtain for a combination of two 
realistic magnetized iron detectors with L = 732 and L = 3000 baselines. However, we 
stress here that the OPERA-like 2 Kton lead-emulsion detector is under construction 
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and should be operational by mid 2006 to exploit the CNGS beam. We can therefore 
study how possible improvements of the detector (from 2006 till the time the Neutrino 
factory will become operational) affect the results of our analysis. In particular we 
remind that a moderate scaling in the lead-emulsion detector mass could indeed be se- 
riously taken into account (due to the rapid increase in the emulsion scanning power). 
The outcome of this study is that an increase in the emulsion mass from 2 to 4 Kton 
eliminates completely the "clone" regions for 13 > 1°. On the other hand, a moderate 
increase in the main background rejection does not seem to improve in a significant 
way the previous results. 

As a final comment we believe that the two-detector types combination should be 
investigated further, even in the case of results comparable to those of a two-baselines 
magnetized iron detector combination, since having detectors with different character- 
istics and systematics can be extremely helpful to ameliorate the 5 measurement. 
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A Perturbative expansion in A6 



If #13 is large enough, we can expand eqs. ([5]) and (Q) in powers of AO for any value of 
A6 in the allowed region of Fig. [l|. 

At first order in Ad, 



(24) 



/(0i3, 13 ) * -A9 fiflra) + 0( A# 2 ) ~ -4A# |g + 0(A# 2 ) , 

g(8i3,0i3) * 1 - AO 13 ) + O(A0*) ~ 1 - Afl + 0(A# 2 ) 

where, for #i 3 e [0°, 13°], / and <? are positive functions of #13. 
For neutrinos we get 

A , L f /- A t L\ ARJ(0 13 ) + cos (5 - ^mL) g(9 l3 ) 

2 ^ V 2 / I S i n ( § — &atmL 

and for antineutrinos we get 

5 _ = _^nL ± U- 5+ A^L. ^ ' ' 7 } . ( 20) 





) + cos 




9(013) 




sin 


f 5 + ^f^) 





Notice that from eq. (|5]), for A6 = 0, we get 5± = 5 or 5± = — 5 ± A atm L. For A6 1 7^ 0, 
we get the two branches of the equiprobability curves. 

If we turn on the AO correction at first order, the results of Fig. ^ are easily under- 
stood with the following argument: consider two different equiprobability curves for 
neutrinos, for two different values of the neutrino energy Ex and E 2 , E 2 > E\. The 
intersections between the two curves are defined by 

S+iEi) = 5+(E 2 ). 



Solving for AO we get two solutions: 

AO = 0, 

AO T A a tm(El) — A a tm(E2) 

^° ~ ^ H(E 2 )-H(E 1 ) ' 

where 

AR + (E)f(0 13 ) + cos (5 + ^mk) 13 ) 



(27) 



H(E) 



sin (s + ^mk} 
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For AO = 0, we get 5 = 5. On the other hand, the second intersection gets an O(A0) 
correction: 

5 = -6 + LA atm (Ei) + A0ff . (28) 

Being H(E) and (for the considered parameters and baselines) the second solution for 
A9 both positive quantities, the second intersection between any two equiprobability 
curves for neutrinos will be displaced towards negative values of 5 and positive values of 
AO. This is precisely what observed in the upper row of Fig. ^|. For the considered set 
of input parameters, baselines and neutrino energy range the same argument applies 
to the antineutrino equiprobability curves also, as it can be seen in the lower row of 
Fig. |. 



B Formulae for r decay and v T N cross-section 



Grouping events in bins of the final muon energy E^, with the size of the energy bin 
depending on the energy resolution AE^ of the considered detector, the number of 
"silver" muon events in the i-th energy bin for the input pair (# 13 , 5) and for a parent 
muon energy E ^ is: 



NU0 13 ,5) = BR(r^fi) 



P±(E V ,9 13 ,6) 



dN^{E^E T ) da MBr) (E r ,E v ) 



dE„ 



dE T 



d$v e (v e )(E u , Efj) 



Ei + AEr, 



dE u 



where ® stands for a convolution integral on the intermediate energy. 
Several comments are in order: 



The differential neutrino flux d^ Ue ^ e )(E v , E^) jdE v has been given in eq. fllPP- 
The v e — > v T oscillation probability has been given in eq. (f2~0|). 
The v T N cross-section 

The (anti)neutrino-nucleon differential cross sections on an isoscalar target, de- 
fined as the average of the differential cross-sections on proton and neutron, can 
be divided in three components: the elastic, the quasi-elastic and the DIS cross- 



section. We refer to 77 for a detailed discussion on the different contributions 



to these components. The DIS cross-section can be expressed in terms of the 
isoscalar structure functions Fi(x) (i = 1, . . . , 5): 



dxdy 



G 2 F ME V 

7i 



y\xy + 



2ME U 



F 1 {x)+ [l-y 



Mxy 
~2E~ 



F 2 (x) 



± 



x v ( 1 - 2 



mt 



y 



AME„ 



F 3 (x) + \xy 



mt 



+ 



4 

mf 



2ME„ AM 2 El 



FAx) 
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where M is the nucleon mass and mi the charged lepton mass; x = Q 2 /2ME u y is 
the DIS Bjorken variable and y — 1 — Ei/E v depends on the fraction of neutrino 
energy carried away by the charged lepton. Notice that the x, y variables defined 
in this Appendix should be interpreted as local variables: outside the Appendix 
the same letters may represent different objects (we are simply following the 
standard notation for the DIS cross-section). 

Making use of the Callan-Gross relation, 2xFi(x) = Fi(x), and of the Albright- 
Jarlskog relations |7^| , F±{x) = 0,xF 5 (x) = F 2 (x), eq. (p9|) depends on two 
independent structure functions, only. Using the quark parton model (see |78| 
and refs. therein, for example) the isoscalar structure functions F 2 {x) and xF 3 (x) 
are: 

x \(u + d + 2s) + (u + d + 2c) 
(u + d + 2c) + (u + d + 2s) 



F?(x) 



FZ(x) 



xF%\ 



x 



xF»(x) 



X 



X 



X 



(30) 



d + 2s) - (u + d 
u + d + 2c) - (u + d- 



2c) 
2s) 



Eventually, the structure functions, given in terms of the parton distribution at 
a definite value of Q 2 [f79| , must be evolved according to the DGLAP equations 
to compute the z/ r iV cross-section for a given neutrino energy. 
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y 



Figure 15: Neutrino and antineutrino - nucleon differential cross section as a function 
ofy for E v = E^. 



In Fig. [T5| we show the dependence of the da/dy on y for a fixed value of E v (in this 



case, En = ES). Notice the strong suppression of the antineutrino cross-section 
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with respect to the neutrino one for high and intermediate r energy (y — ► 0). 
The integration limits for the (x, y) variables are reported in Fig. [16] for different 
neutrino energies [3^ • 




o.i 04 o.j at. Q7 as. «.t 



Figure 16: Kinematical bounds in the x, y plane for different neutrino energies, j^Qj. 



The r decay rate 

The muon distribution is given, in the r rest frame, by the following expression: 



d 2 N fl ± 
dx'd cos 9 



[f {x',R)TV T h{x\ R) cos#] 



(31) 



where R = m^/m T , V T is the average r polarization along the r direction in the 
laboratory frame, x' = 2E u /m T (not to be confused with the Bjorken x variable 
defined for the cross-section) and 9 is the angle between the muon momentum 
vector and the r spin direction. The two /-functions are: 



fo(x',R) = 2Vx' 2 - 4R 2 [-4R 2 + 3(1 + R 2 )x - 2x' 2 \ (32) 
f\{x',R) = 2{x' 2 -AR 2 ) [l + 3R 2 -2x'] (33) 

and they reduce to the standard fo(x') and fi(x') functions for the and e 



flux in the R — > limit |81| . The boost in the laboratory frame is given by the 



31 



following relations: 

nJ ft . 

x' 2 - 4R 2 



x> 8 

V = 

y 2 2 



1/2 

cos 6 



Px' + (x' 2 -AR 2 ) 1/2 cos9 / o a\ 

COS 6' = -pr (34) 

[{X' + [3(x> 2 - 4i? 2 ) 1 /2 cos Q) _ 4(1 _ £2)^2] 1/2 

where y' = E^/Er and 9' is the angle boosted in the lab frame. 
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